Intermode Dephasing in a Superconducting Stripline Resonator 
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We study superconducting stripline resonator (SSR) made of Niobium, which is integrated with a 
superconducting interference device (SQUID). The large nonlinear inductance of the SQUID gives 
rise to strong Kerr nonlinearity in the response of the SSR, which in turn results in strong coupling 
between different modes of the SSR. We experimentally demonstrate that such intermode coupling 
gives rise to dephasing of microwave photons. The dephasing rate depends periodically on the 
external magnetic flux applied to the SQUID, where the largest rate is obtained at half integer 
values (in units of the flux quantum). To account for our result we compare our findings with 
theory and find good agreement. 



INTRODUCTION 

Dephasing is the suppression process of quantum co- 
herent effects due to coupling between a quantum system 
and its external environment [1]. A Kerr nonlinearity in 
an electromagnetic resonator may lead to dispersive in- 
termode coupling, which in turn may result in dephasing 
of photons [2, 3]. Such a coupling mechanism can also 
be exploited to allow quantum non demolition (QND) 
detection of single photons [2, 4-7]. A Kerr nonlinear- 
ity exists in superconducting stripline resonators (SSR) 
due to the effect of kinetic inductance. However, the re- 
sultant intermode coupling is typically far too weak to 
allow any significant dephasing [8]. On the other hand, a 
much stronger Kerr nonlinearity can be achieved by in- 
tegrating a superconducting interference device (SQUID) 
with the SSR [9-11]. External magnetic flux can be em- 
ployed in these devices to modulate both the linear and 
nonlinear contributions to the inductance of the SQUID, 
which in turn allows external control of both the reso- 
nance frequencies and the strength of Kerr nonlinearity 
respectively. The enhanced Kerr nonlinearity also pro- 
vides strong coupling between different modes in the res- 
onator that causes dephasing of one mode (called the 
system mode) when another one (the detector mode) is 
externally driven at relatively high powers. 

Here, we employ such a configuration consisting of a 
Niobium SSR and incorporating a SQUID device hav- 
ing a nanobridge in each of its two arms. We moni- 
tor the resonance lineshape of one of the modes of the 
resonator (the system mode) as we simultaneously drive 
another one (the detector mode). We find that a signif- 
icant broadening of the resonance lineshape of the sys- 
tem mode occurs in the same region where the response 
of the detector mode, which is measured simultaneously, 
becomes strongly nonlinear. We provide theoretical evi- 
dence to substantiate our hypothesis that the underlying 
mechanism responsible for the observed broadening is in- 
termode dephasing. The ability to externally control the 
strength of the intermode coupling, which in turn con- 



trols the dephasing rate, makes our device an ideal tool 
for experimentally studying fundamental issues related 
to the quantum - classical transition [1] . 
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FIG. 1: (a) The device and the experimental setup. The 
SSR is made of two identical stripline sections of length 
Zt/2 = 110 mm, each having inductance Lt and capaci- 
ta nce Ct per unit length and characteristic impedance Zt = 
•y/ Lt I Ct = 50 Q. The stripline sections are connected by 
a SQUID consisting of nanobridge type weak links. All mea- 
surements are done at liquid Helium temperature 4.2 K, where 
the device is placed inside a copper package, which is inter- 
nally coated with Nb (to reduce surface resistance and to pro- 
vide magnetic shielding) . (b) SEM micrograph of the SQUID 
(tilted view). The loop area is 39x39/im 2 . (c) SEM micro- 
graph of the nanobridge. 
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Fig. (1) schematically shows the device. The SSR [11- 
13] comprises two identical stripline sections connected 
by a SQUID. A nanobridge [14-17] on each arm of the 
SQUID loop serves as a weak link. The critical currents 
of the nanobridges are denoted by I c \ and I C 2 respec- 
tively. Both nanobridges are assumed to have the same 
capacitance Cj. The self inductance of the loop is de- 
noted by A. A feedline, which is weakly coupled to the 
SSR, is employed to deliver the input and output mi- 
crowave signals. The experimental setup is presented in 
subplot (a) of Fig. (1). 

The fabrication process starts with a high resistivity 
Si substrate coated with SiN layers of thickness 100 nm 
on both sides. A 150 nm thick Nb layer is deposited 
on the wafer using magnetron DC sputtering. Then, e- 
beam lithography and a subsequent liftoff process are em- 
ployed to pattern an Al mask, which defines the SSR and 
the SQUID leads. The device is then etched using elec- 
tron cyclotron resonance system with CF4 plasma. The 
nanobridges are fabricated using FEI Strata 400 Focus 
Ion Beam (FIB) system [18-22] at accelerating voltage 
of 30 kV and Ga ions current of 9.7 pA. The outer di- 
mensions of the bridges are about 150 x 50 nm. However, 
the actual dimensions of the weak-links are smaller, since 
the bombarding Ga ions penetrate into the Nb layer, and 
consequently, suppress superconductivity over a depth es- 
timated between 30 nm to 50 nm [22, 23]. 

EFFECTIVE HAMILTONIAN 

The effective Hamiltonian of the closed system consist- 
ing of the SSR and the SQUID , expressed in terms of 
the annihilation and creation operators Ax, A[, A3 and 
A\ of the first and third modes respectively, is given by: 

H eS = nwxNx + fmj 3 N 3 + V iu 

+ HKxNl + h\x, 3 NxN 3 . (1) 

where Nx — A\Ai and N 3 = A' 3 A 3 are number operators, 
V in = fr\/2~7nfei n [er^Ax + e^A^j represents the ex- 
ternal driving, 7^ is the coupling constant between the 
first mode and the feedline, b™ is the amplitude of the 
driving pump tone which is injected into the feedline to 
excite the first mode, and where u> p is its angular fre- 
quency. Full Derivation of the Hamiltonian is given in 
appendix A. The last two terms represent the Kerr non- 
linearity term of the first, externally driven (detector) 
mode and the intcrmodc coupling between the first and 
the third (system) modes, respectively. The coefficients 
u>x, <^3, Kx and Aj.,3, which are calculated in appendix A., 
depend periodically on the external flux <& x with period 
$o- The flux dependence of lo\ and W3 can be attributed 
to the inductance of the SQUID, which is proportional 
to the second derivative of £q with respect to /, where Sq 



is the ground state energy of the SQUID. On the other 
hand, both the Kerr nonlinearity K\ and intermode cou- 
pling A1.3 coefficients are proportional to the nonlinear 
inductance of the SQUID [24], which in turn is propor- 
tional to the fourth derivative of Sq with respect to /. 



RESONANCE FREQUENCY SHIFT 



Figure (2) shows measurements of the reflection coeffi- 
cient I Sii| (Sn is the ratio between the reflected outgoing 
and the injected incoming amplitudes in the feedline) of 
the first 3 modes of the resonator as a function of fre- 
quency and externally applied flux $ a . The sketches on 
the right hand side show the current waveform of each 
mode. For the first and the third modes, Sxx is found 
to be a periodic function of <f> x with period where 
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FIG. 2: (Color online) Reflection coefficient |Sn| vs. fre- 
quency and external flux for the first 3 modes of the SSR. 
A change of 4.8 mA in the external current corresponds to a 
change of $o in the magnetic flux. 



3 



$o = h/2e is the flux quantum. On the other hand, 
the 2nd mode, which is decoupled from the SQUID since 
its current waveform has a node at the location of the 
SQUID, does not exhibit a flux dependence. Note that 
the data in Fig. (2) is obtained by sweeping the mag- 
netic flux upwards. However, as can be seen from Fig. 
(3a), in which the resonance frequency f± — lui/2tt of the 
first mode is measured versus both increasing (blue) and 
decreasing (red) magnetic flux, the response is hysteretic. 

The solid black line in Fig. (3a) is obtained by nu- 
merically evaluating the resonance frequency fi = lui/2tt 
using Eq. (7) in appendix A. For the parameters that 
are used in the calculation for this case (see figure cap- 
tion), the SQUID can be either monostable or bistable 
depending on Q x . Consequently, sharp transitions oc- 
cur near the values of Q x corresponding to a bound- 
ary between these regions, and as a result, the response 
is hysteretic. Interestingly, as the input power is in- 
creased the response becomes non-hysteretic, as can be 
seen from Fig. (3b), which shows a measurement of /i 
at Pj n = — 71 dBm. Theoretically, this behavior is ac- 
counted for by assuming that the value of the screening 
parameter /3l = 27rA/ c /$ , which is proportional to the 
average critical current I c = (I c i + I C 2) /2, is significantly 
lower for this case (0.15 instead of the value 7.4, which 
was used to fit the data for p n = —96 dBm). To account 
for this behavior we discuss in appendix B. the possibility 
that local heating of the nanobridges is responsible for the 
drop in I c at elevated input powers. Assuming that the 
heat is mainly dissipated down into the substrate rather 
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FIG. 3: (Color online) The resonance frequency shift A/i of 
the first (detector) mode vs. applied flux for two different 
values of P n . The flux is first swept upwards (blue line) and 
than downwards (red line) . The black solid lines represent the 
theoretical calculation of A/i using the following parameters: 
(3 L = 7.4 for P in = -96 dBm, /3 L = 0.15 for P in = -71 dBm, 
and P1/P2 = 3 for both cases. 



than along the film, we estimate that the temperature 
rise for P n = —70 dBm is 4K. This rough estimation in- 
dicates that heating may indeed play an important role, 
and may be held responsible for the apparent drop in the 
critical current. 



INTERMODULATION AND INTERMODE 
DEPHASING 

Both nonlinear terms in the Hamiltonian H. c g (1) play 
an important role as p n is increased. The effect of the 
Kerr nonlinearity can be sensitively observed by employ- 
ing intermodulation (IM) characterization [24]. In this 
method, in addition to the relatively strong pump tone 
at frequency u> p , which is used to drive the first mode to 
any desirable operating point, another tone, called sig- 
nal, which has a much smaller power P s .i n and a nearby 
frequency uj p + Slu (Slo is much smaller than the res- 
onance width), is also injected simultaneously into the 
fccdline. Due to Kerr nonlinearity these two inputs may 
mix in the resonator and produce tones of IM products. 
Typically, the largest IM products are the output sig- 
nal at frequency oj p + 5lo and the output idler at fre- 
quency cl> p — 5lo. The two corresponding gain factors, 
namely the signal gain G s = P s ,out/P,in and the idler 
gain G; = Pi, ou t/P s ,in, where P s , out and P iiOUt are the 
powers of the output signal and output idler tones re- 
spectively, were evaluated in Ref . [24] . Panel (b) of Fig. 
(4) presents a color map showing IM characterization of 
the first (detector) mode, which was obtained using a 
spectrum analyzer. The powers of the injected pump 
and signal tones in the IM measurement are —62.1 dBm 
and —81 dBm respectively. Both G s and G{ periodically 




I Snl [dB] 



-100 -90 -80 -70 -60 -50 

< b > Reflected Power [dBm] 



FIG. 4: (Color online) IM characterization of the detector 
mode [panel (b)] and |Sn| measurements of the system mode 
[panel (a)]. Largest idler gain as well as highest dephasing 
rate is obtained at half integer values of the externally applied 
flux. 
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FIG. 5: (Color online) Signal gain G s , idler gain Gi, and nor- 
malized dephasing rate T v . Experimental results are shown 
in panels (al), (bl) and (cl), whereas theoretical predictions 
are shown in panels (a2), (b2) and (c2). The following device 
parameters were taken in order to evaluate G s and Gi theo- 
retically using Eqs. (82) and (83) of Ref. [24] [panels (a2) 
and (b2) respectively)], and to evaluate T v using Eq. (70) of 
Ref. [8] [panel (c2)]: ( 7f3 + 7d3 ) /lui = 5000, 7 f3/7d3 = 0.15, 
id = 1.5 and I C 2 = 4.5 ^iA. Note that the effect of non- 
linear damping [24] is disregarded. 



oscillate as a function of the external current. This be- 
havior is seen more clearly in panels (al) and (bl) of Fig. 
(5) which exhibit G s and Gi versus $ K /$o- 

The intermode coupling term in the Hamiltonian (1) 
can be exploited to continuously measure the number of 
photons in the system mode by externally driving the de- 
tector mode [2, 3]. Such a measurement scheme is char- 
acterized by the time it takes to resolve adjacent number 
states of the system mode. Significant dephasing occurs 
when this time scale is made comparable or shorter than 
the lifetime of photons in the system mode. Theoret- 
ically, dephasing of photons in the system mode is ex- 
pected to give rise to a resonance frequency shift and 
to broadening of the resonance line shape of the power 
reflection coefficient |5n (w)| 2 , which is given by [25] 



|5nH| 2 = l- 



7f37d3 



4 7 t, 



7f3 + 7d3 7t 2 ot + (w - ua) 



(2) 



where (D3 is the shifted angular resonance frequency, the 
total width is given by 7t ot = 7f3 + 7d3 + 1/%,, where 7f3 
denotes the coupling constant between the system mode 
and the feedline, and 7 d3 denotes the damping rate of the 
system mode and l/r^ is the dephasing rate of photons 
in the system mode. 

Simultaneously with the IM characterization, we also 
measure the resonance line shape of the third mode using 
a network analyzer. A very low input power of —101 
dBm is employed to avoid any nonlinear response of the 



third mode. As can be seen from the results, which are 
presented in panel (a) of Fig. (4), the measured reflection 
coefficient \Sn (w)| periodically oscillates as a function 
of Q x . Fitting the experimental data to Eq. (2) yields 
the normalized dephasing rate r v = 1/ ( 7 f3 + 7H3) t v . As 
can be seen from Fig. (5), at the same points where Gi 
peaks [panel (bl)], namely for half integer values of the 
external flux, a strong peak is found in T v [panel (cl)]. 
At these points, the value of T v exceeds unity, namely, 
the dephasing rate becomes larger than the system mode 
decay rate. 

To account for the experimental results, we employ 
Eqs. (82) and (83) of Ref. [24] to calculate the gain 
factors G s and G; respectively, and Eq. (70) of Ref. [8] to 
calculate the normalized dephasing rate T v . The results, 
given in panels (a2), (b2) and (c2) of Fig. (5) yield fairly 
good agreement with the experimental data [panels (al), 
(bl) and (cl)]. The device parameters that were used 
in the calculation are listed in the figure caption. The 
flux dependence of the gain factors G s and G; and that 
of the normalized dephasing rate T v can be attributed 
to the periodic flux dependence of the parameters oji, 
U3, K\ and Ai,3 of the Hamiltonian (1). Both nonlinear 
parameters K\ and A13 peak at half integer values of 
the external flux. Consequently, both Gi, which can be 
considered as a measure of the strength of nonlinearity, 
and T v , which strongly depends on A1.3, obtain their 
largest values at these points. 



CONCLUSION 

Integrating a SQUID having large nonlinear induc- 
tance with an SSR leads to strong IM distortion and 
strong intermode coupling. In the present paper we have 
exploited these effects to study a novel mechanism of de- 
phasing of microwave photons that can be externally con- 
trolled. The same intermode coupling that is responsible 
for the observed photon dephasing can also be exploited 
for single photon detection[2, 3]. In future experiments 
several improvements, such as increasing the nonlinear 
coupling, as well as reducing the temperature and using 
lower noise pre-amplifier should allow detection of single 
microwave photons. 
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APPENDIX A: DETAILED DERIVATION OF 
THE EFFECTIVE HAMILTONIAN 

The effective Hamiltonian of the closed system com- 
prising the SSR and the SQUID [12, 26] is found using 
the same method that was previously employed in Refs. 
[13, 26]. Here however, we relax the assumption that the 
self inductance of the SQUID loop is small, and also the 
assumption that both junctions have the same critical 
currents. On the other hand, we assume that the induc- 
tance of the SQUID, which is denoted as L$ 7 is much 
smaller than the total inductance of the stripline Lt^t- 
This assumption can be justified by considering the fact 
that the measured angular resonance frequencies lj„ of 
the first 3 modes (n e {1, 2, 3}) for all values of & x (see 
Figs. 2 and 3 in the paper body) are very close to the 
values expected from a uniform resonator having length 
It, namely nwx, where ivt — tt/ItVLtCt- Moreover, 
the normalized flux- induced shift Au) n /nujT in the angu- 
lar resonance frequency of the first 3 modes is quite small 
and never exceeds 10 -3 . Both observations indicate that 
the ratio Ls/LtIt can indeed be considered as a small 
parameter. 

The resultant Hamiltonian of the closed system is given 
by H — Hssr + Us {I), where Hssr is the SSR Hamil- 
tonian and where Hs (I) is the SQUID Hamiltonian, 
which depends on the current / at the center of the SSR, 
namely, the current flowing through the SQUID. In terms 
of annihilation (A± and A 3 ) and creation (A\ and A\) 
operators for the first and third modes of the SSR re- 
spectively the Hamiltonian Hssr can be expressed as 

Hssr = fiur (JVi + 3iV 3 ) + M„ , (3) 
where Ni — A\A X and N 3 = A^ are number operators, 



-iuj-pt 



pt Ai + e 1 



(4) 



represents the external driving, -fn is the coupling con- 
stant between the 1st mode and the feedline, 6 1 " is the 
amplitude of the driving pump tone, which is injected 
into the feedline to excite the first mode, and uj p is its 
angular frequency 



The kinetic inductance of the nanobridges 

The Hamiltonian for the SQUID depends on the prop- 
erties of the nanobridges. Due to the Ga ions im- 
planted in the outer layer of the Niobium during the 
FIB process and the consequent suppression of super- 
conductivity in that layer [22, 23], the weak links are 
treated as variable thickness nanobridges. The behav- 
ior of such a nanobridge is strongly dependent on the 
ratio l/£ [14, 17, 27-32], where I is the bridge length 
and £ is the coherence length of the Cooper pairs. The 



coherence length £ depends also on the temperature of 
the bridge. In the dirty limit £ is given by £(T) = 



0.852yJt lf(T c /T-l) 1 [14], where £ is the size of 
the cooper pair and If is the mean free path [33, 34]. 
The current-phase relation (CPR) of the bridges is pe- 
riodic with respect to the gauge invariant phase 6 across 
the bridge. When l/£(T) <C 1, the nanobridge be- 
haves like a regular Josephson junction (JJ) with a si- 
nusoidal CPR[35]. However, as the ratio l/£(T) becomes 
larger, the CPR deviates from the sinosoidal form and 
can also become multivalued [14]. In case the CPR is 
not multivalued the bridge can be approximately con- 
sidered as a JJ having an extra kinetic inductance Lk 
. The effect of the kinetic inductance can be taken into 
account by replacing the screening parameter of the loop 
0l = 27rA/ c /$o by an effective one given by (3l + A/3, 
where A/3 = 2ttLkI c /$o- 

In order to estimate A/3 we use Eqs. (47)-(49) and the 
data in Fig. 5ofRef. [15]. For?/£ = 1.7 the bridges' con- 
tribution is A/3 ~ 1 . As we will discuss below, both /3l 
and A/3 depend on the injected power P m that is used 
to excite the resonator due to a heating effect. How- 
ever, for all values of P m that were used in our experi- 
ment, we estimate that the ratio A/3//3l never exceeds 
the value 0.5 and thus the effect of kinetic inductance can 
be considered as small. Furthermore, the CPR remains 
a single valued function in the entire range of parame- 
ters that is explored in our experiments. Consequently, 
the nanobridges can be treated as regular JJs to a good 
approximation. 

The SQUID Hamiltonian 

In the following derivation we treat the nanobridges 
as regular JJs. We consider the case where the criti- 
cal currents of both nanobridges are J c i = J c (1 + a) and 
I C 2 = I c (1 — a ) respectively, where the dimensionless pa- 
rameter a characterizes the asymmetry in the SQUID. 
The Hamiltonian for the SQUID, which is expressed in 
terms of the two gauge invariant phases <5i and <5 2 across 
both junctions, and their canonical conjugates p\ and P2, 
is given by 



n s (i) 



En 



+ E u{5 u S 2 ;I) , (5) 



where u p i = \Jl c /Cj<&n is the plasma frequency, E — 
<f>nI c /-K is the Josephson energy, and the dimensionless 
potential u is given by [36] 



( &\ — $2 i rti I 

1 + a) cos<5i + (1 — a) cos(5 2 \ 2 ' # Q J 



(61+62)! C(6i+6 2 f 



Air 
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(6) 
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where ( = ^ /^IcL t It- 

Adiabatic approximation 

Due to the extremely small capacitance Cj of both 
nanobridges [37] , the plasma frequency cu p i of the SQUID 
is estimated to exceed 1 THz. Thus, the effect of the 
SQUID on the SSR, which has a much slower dynam- 
ics, can be treated using the adiabatic approximation 
[6, 38]. Formally, treating the current J as a parame- 
ter (rather than a degree of freedom), the Hamiltonian 
Hs can be diagonalized Hs \k (I)) = Sk (I) \k(I)), where 
k = 0, 1, 2, and (k (I) \l (I)) = 5 M . To lowest order in 
the adiabatic expansion the effective Hamiltonian gov- 
erning the dynamics of the slow degrees of freedom cor- 
responding to the fast part of the system occupying the 
state \k{I)} is given by H% = Hssk + e k {I) [39, 40]. 
Furthermore, in the limit where the thermal energy k^T 
is much smaller than the typical energy spacing between 
different levels of Hi (~ hu> p \) one can assume that the 
SQUID remains in its current dependent ground state 
|0 (/)) . For most cases this assumption is valid for our ex- 
perimental parameters. It is important, however, to note 
that when the externally applied magnetic flux is close 
to a half-integer value (in units of $o), namely, when 
<& x ~ (n + 1/2) $0, where n is integer, this approxima- 
tion may break down. Near these points the potential 
u may have two different neighboring wells having sim- 
ilar depth. Consequently, near these points, the energy 
gap between the ground state and the first excited state 
can become much smaller than htOpi. On the other hand, 
the ratio between the height of the barrier separating 
the two wells (~ E ) and the energy spacing between 
intra-well states (~ Huj p \) is typically E /huj p \ ~ 100 for 
our samples. Since the coupling between states localized 
in different wells depends exponentially on this ratio, we 
conclude that to a good approximation the inter- well cou- 
pling can be neglected. Moreover, in the same limit where 
Eo/hujpi 3> 1, one can approximate the ground state en- 
ergy £o by the value of Equ at the bottom of the well 
where the system is localized. 

The current I at the center of the SSR can readably 
be expressed in terms of the annihilation and creation 
operators A\, A\ A$ and A\. This allows expanding the 
current dependent ground state energy £o (/) as a power 
series of these operators. In the rotating wave approx- 
imation oscillating terms in such an expansion are ne- 
glected since their effect on the dynamics for a time scale 
much longer than a typical oscillation period is negligibly 
small. Moreover, constant terms in the Hamiltonian are 
disregarded since they only give rise to a global phase 
factor. In the present experiment the 1st SSR mode is 
externally driven, and we focus on the resultant dephas- 
ing induced on the 3rd mode. To that end we include in 
the effective Hamiltonian of the closed system in addi- 



tion to the linear terms corresponding to the 1st and 3rd 
modes, also the Kerr nonlinearity term of the 1st mode, 
which is externally driven, and also the term representing 
intermode coupling between the 1st and the 3rd modes 
[seeEq. (1)]. 

The angular resonance frequency shift of the 1st and 
the 3rd modes, which is given by 

- m T _ lo 3 - 3a> T _ d 2 {sq/Eq) 
u T ~ 3uj t ~ d{I/I c f ' 

can be attributed to the inductance of the SQUID, which 
is proportional to the second derivative of £o with respect 
to /. On the other hand, the Kerr nonlinearity, which is 
given by 

Ki C 2 ^i d A (so/E ) 

2E Q d(I/I c f ' 1 ' 

and the intermode coupling, which is given by A13 = 
9Ki , can both be attributed to the nonlinear inductance 
of the SQUID [24], which is proportional to the fourth 
derivative of £0 with respect to /. 

Evaluation of ui\, uz, K\ and A13 in the limit /3l <S 1 

The evaluation of the parameters wi, u>s, K\ and A13 
generally requires a numerical calculation. However, an 
analytical approximation can be employed when <C 1. 
In this limit the phase difference 82 — Si is strongly con- 
fined near the value 2tt& x /&o, as can be seen from Eq. 
(6). This fact can be exploited to further simplify the 
dynamics by applying another adiabatic approximation, 
in which the phase difference S2 — Si is treated as a 'fast' 
variable and the phase average 5+ = (Si + S 2 ) /2 as a 
'slow' one. To lowest order in the adiabatic expansion 
one finds that for low frequencies u) -C to p \, namely in 



Forward Flux Sweep Backward Flux Sweep 
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i a \ Frequency [GHz] ^ Frequency [GHz] 

FIG. 6: Measured |6ii| at input power Pj n = — 95dBm for 
forward (a) and backward (b) magnetic flux sweep. In this 
sample (3l = 20, and the response is highly hysteretic. 
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the region where the impedance associated with the ca- 
pacitance of the JJs is much larger in absolute value in 
comparison with the impedance associated with the in- 
ductance, the SQUID behaves as a single JJ having crit- 
ical current given by [41] 



I s = 2W1 - (1 - a 2 ) sin^ (tt^/Sq) 



(9) 



Note that this approximation may break down when 
Q x ~ (n + l/2)<I>o unless the asymmetry parameter a 
is sufficiently large. The relatively large value of a in 
our device (a ~ 0.5) ensures the validity of the above 
approximation. Using this result, it is straightforward to 
obtain the following analytical approximations: 



d 2 (eo/E ) 
d(I/Ic) 2 
d 4 (eo/Eo) 



irl s 
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(10a) 
(10b) 



d {I I let 

which can be used to evaluate all the terms in Eq. (1) 



APPENDIX B: HYSTERETIC RESPONSE AND 
HEATING OF THE NANOBRIDGES 

As we discuss in the paper, the resonator exhibits hys- 
teretic response to magnetic flux when the input power 
is relatively low. Such a behavior occurs, as can be 
seen from Eq. (6) above, when the screening parame- 
ter (3l is sufficiently large to give rise to metastability 
in the dimensionless potential u. A fitting of the model 
to the experimental data shown in Fig. 3(a) of the pa- 
per yields a value of (3l = 7.4. Another example of hys- 
teretic response is shown in Fig. 6 below that shows data 
taken with another sample, which was fabricated using 
the same process that is described in the first section. 
The larger critical current in that sample yields a larger 
value of the screening parameter (5l = 20. 

As is mentioned in the paper, as the input power is in- 
creased the response becomes non-hysteretic. The grad- 
ual transition between the hysteretic region to the non- 
hysteretic one is seen in Fig. 7 below, which shows the 
difference in the measured resonance frequency of the 
first mode obtained from increased flux sweep (fu n c) and 
decreased flux sweep (/idee) at different input powers. 
Dark blue in the color map corresponds to no difference, 
namely to monostable regions, whereas in the red re- 
gions, where a large difference is observed, the system 
is bistable. As can be clearly seen from the figure, the 
bistable regions shrink as the input power is increased. 
The experimental results suggest that the critical current 
of the nanobridges drops as the input power is increased, 
and consequently the response becomes non-hysteretic 
due to the resultant smaller value of the screening pa- 
rameter Pi,. We hypothesize that the drop in the critical 



current occurs due to heating of the nanobridges by the 
input power. 

To estimate the effect of heating, we assume the case 
where the substrate is isothermal and that the heat is 
mainly dissipated down into the substrate rather than 
along the film [42]. Moreover, we assume that most of 
the externally injected power into the resonator is dissi- 
pated near the nanobridges, where, the current density 
obtains its largest value. By estimating the heat transfer 
coefficient per unit area between each nanobridge and 
the substrate beneath it (100 nm SiN on top of high- 
resistivity Si) to be k ~ lWcm _2 K _1 [43, 44] and the 
area of the nanobridge to be A ~ (50nm) 2 one finds that 
the expected temperature rise for P m = —70 dBm is 
AT = Pi n /AK ~ 4 K. 

Since heating is produced by AC current flowing 
through the nanobridges, it is important to estimate 
also the thermal rate, which characterizes the inverse of 
the typical time scale of thermalization, and is given by 
7t = Ak/C, where the heat capacity C of the nanobridge 
is given by C = C v Ad, C v is the heat capacity per unit 
volume, and d is the thickness of the superconducting 
film. Using the estimate C v ~ 10 -3 Jcm~ 3 K _1 [44] one 
finds 7t — 0.1 GHz. Since the frequency of the AC heat- 
ing current is 1-2 orders of magnitude higher, we con- 
clude that to a good approximation the temperature of 
the nanobridges can be considered as stationary in the 
steady state. 




l f „„ c - f ,dJ[ GHz l 



FIG. 7: The difference between the measured resonance fre- 
quencies obtained in the increasing flux sweep (/line) and the 
decreasing flux sweep (/idee) of the first (detector) mode. The 
dark blue areas correspond to monostable regions, namely, the 
same resonance frequency is measured for both the increased 
and decreased sweep. The red indicates the regions where the 
system is bistable. 
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